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Abstract 


The dynamics of the motion of the Small Expendable Deployer System (SEDS) is studied by using 
a simplified model in which no external forces (except the gravity gradient field) are applied on the 
tethered body and the tether is assumed massless. The dynamics of SEDS operation is modeled as 
a sequence of two phases: the deployment phase and the swing phase. For the first one the velocity 
dependent forces are found to force the tether forward from the local vertical. When the deployment 
ends, Coriolis effects vanish and the swing phase begins, which is characterized by a wide free libration. 
The time duration as well as velocity, acceleration and tension of the tethered body are estimated for 
both deployment and swing phases. 
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Introduction 


This report describes the work done at NASA Langely Research Center during the period from 
August 4, 1988, to September 2, 1988, to model and study the dynamics associated with the use 
of tethered satellites deployed downward from a Space Transportation System (STS) or an ELV 
(Expendable Launch Vehicle) in order to obtain aerothermodynamic and atmospheric science data, 
as well as data about the dynamics of the tethered system itself. 
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The “Tethered Wind Tunnel” concept has been established and has been shown to be a valuable 
adjunct to other proposed flight programs to obtain aerothermodynamic and atmospheric data at 
altitudes between 50 and 150 km [1]. The second Tethered Satellite System (TSS-2) [2], the Shuttle 
Tethered Aerothermodynamic Research Facility (STARFAC) [3] and the Shuttle Continuous Open 
Wind Tunnel (SCOWT) [4], in particular, have the additional capability of obtaining steady-state 
data over the period of one or more orbits at a sustained altitude. However, both TSS/SCOWT and 
STARFAC are relatively large and sophisticated systems. In addition, the design of the satellite and 
of the deployer systems and the mission profiles are dictated by safety and operational constraints of 
STS, so that flight opportunities may be limited. This circumstance led to the definition of a simplified 
system consisting of a sacrificial satellite to be deployed as a secondary payload from an ELV (such as 
a Delta II), or as a primary payload on a smaller ELV. This simplified system is generally referred to as 
the Tethered Dynamics Explorer (TDE), or Tethered Atmospheric Probe (TAP) [5], a non-recoverable 
smaller satellite for obtaining data describing deployment dynamics, tether deployment rate and shape, 
tether tension, acceleration of the satellite, and other engineering and science data. These data might 
also be useful in defining the TSS-2 mission; however, the flight profiles will have little similarity to 
that of the TSS, hence, they will be applicable only in a generic sense. 

This work, which was carried out in a cooperative NASA/ASI (Italian Space Agency) program 
about deployment and retrieval of tethered satellites, dealt with the study of the dynamics of the 
above described systems. It moved in two directions. On one hand, a detailed modeling, performed by 
using sophisticated general-purpose computer codes, has been approached. In particular, the use of the 
SKYHOOK program [6] has been studied. On the other hand, in order to define the order of magnitude 
of the main parameters affecting the dynamics of the DELTA II — Small Expendable Deployer System 
(SEDS) Flight Experiment [7] , a simplified form of the motion equations has been derived and studied. 
Results and recommendations are presented. 

In this report, part I describes the study performed using the SKYHOOK program; part II shows 
the derivation of the equations used to analyze the SEDS dynamics; and part III summarizes the results 
and includes recommendations. 

Part I - Use of the SKYHOOK Computer Code 

The study of the detailed dynamics of tethered systems dealt with a more in-depth knowledge of 
the SKYHOOK program [6], originally developed by the Smithsonian Astronautical Observatory, with 
the view of modifying codes to run on a 80386 class 20 MHz desktop computer. The SKYHOOK is 
a program of great generality which has the capability of analyzing a broad range of tether related 
scenarios. The STS and the satellite are represented by mass points, while the tether may be 
included in either mass or may be represented separately as a series of independently spaced mass 
points. Both thermal and elastic expansion of the tether, as well as solar and lunar perturbations, 
electrodynamic forces acting on the tether, and the effects of gravitational fields may be included. 
Moreover, atmospheric drag on the tether and the oblateness of the Earth may also be incorporated into 
the calculation. Analytical Mechanics Associates Inc. has modified the original version of SKYHOOK 
to include lower altitude calculations and a deployment routine specific to the STARFAC which more 
accurately deploys and retains the satellite at a target altitude for the required number of orbits. 

The SKYHOOK mainly consists of the following routines: 

SKYHOOK - Main driver program which calls SKYIO to read data describing the run, and SKYIN2 
which controls the integration of the equations of motion of the tether system. 

CONTRL - Computes the control-law tension on the wire segment being deployed. It is called by 
TENSION. 

DIFSUB - Integrates a set of N ordinary differential equation of first order, in our case the equation 
of motion plus temperature and charge equations. This routine, originally written 
by Gear, uses either an Adams-Moulton predictor — corrector method or a multistep 
method suitable for stiff systems. DIFSUB is called by SKYTN2. 
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- Compute the right hand side of equations to be integrated. In particular, DIFFUN 
computes the sum of all forces acting on each segment of tether, the wire temperature 
(sub WIRTEMP) and the charge rate for electrodynamics (sub CHGRATE). DIFFUN 
is called by SKYIN2 and DIFSUB. 

- Points each data frame at requested frequency, specified by the variable DELT. It is 
called by SKYIN2. 

- Controls the computation of the initial conditions for a mass about to be deployed. It 
is called by SKYHOOK. 

- Set up the initial state vector for a mass about to be deployed. INIT2 is called during 
the deployment phase by INITAL. 

- Read (and print) mass properties. It is called by SKYIO. 

- Read (and print) the harmonic coefficients for the Earth gravity. READCS is called by 
SKYIO. 

- This subroutine is the “true” SKYHOOK iptegrator. It controls the integration of the 
equations by calling DIFSUB, DIFFUN, DMPR. SKYIN2 contains two main loops. 
The inner loop performs the integration of the equations using time step H, the outer 
loop, ending every DELT time step, evaluates the different quantities to be printed 
and recomputes the nominal length of the wire segment just deployed according to the 
Hook law (sub NEWLO). SKYIN2 is called by SKYHOOK. 

- Controls the reading of all input parameters at the beginning of the run. It is called 
by SKYHOOK. 

- Controls the computation of the tension on a section of the wire. If the tether is being 
deployed or retrieved, TENSION calls the proper control law tension for the segment 
being reeled out or pulled in, using CONTRL or RETRV, respectively. Tensions on the 
other sections are computed by using Hook law and damping with subroutine TETHER. 
TENSION is called by DIFFUN. 


In order to better understand the overall operating conditions of the SKYHOOK program, and in 
particular to handle input, output and integrator modules, two runs have been carried out and the 
relative outputs studied. In the first run the downward deployment of a 20 km tether front an orbiter 
at 220 km of altitude has been simulated, using a 3-mass tether model. In the second run, a 4-mass 
tether model has been used to simulate the tether deployment from an orbiter at the same altitude of 
220 km, down to 90 km. In both runs an equatorial orbit has been assumed. 


Part II - Dynamics of Small Expendable Deployer System (SEDS) 


The SEDS is based on two main concepts [8]: low tension deployment (which causes a large libration) 
and use of an expendable tether and satellite. In a typical SEDS operation, the payload is initially 
ejected from the STS by using a spring. Then the payload drifts away and pulls the tether out under 
very low tension. Actually, the system needs small tension adjustments to maintain the deployment 
schedule and a final braking phase reducing the deployment velocity at the end of the 20 km tether to 
prevent rebound of the end mass. At the end of the deployment phase a large in-plane (pitch) angle 
is obtained and, the effects related to the deployment rate vanish, a wide free libration ensues (swing 
phase). Finally, when the payload is about at a vertical position, the tether is cut and the payload is 
released. A detailed discussion about the main advantages related to the SEDS concept (low-tension 
deployment and expendable tether) is out of the aim of this report. It may be found in ref. [8]. In 
the following, in order to evaluate the order of magnitude of the main parameters characterizing the 
dynamics of SEDS deployment and swing phase, a rather simplified form of the equations of motion 
will be derived. It is believed that many of the most important aspects of the behavior of the SEDS 
dynamics can be predicted without resorting to complicated modeling and sophisticated computer 
codes. 


4 



II.l - DEPLOYMENT PHASE 


For this study, the following equations of motion [9] are dcrivated: 

=3 cos 2 9 cos 2 ?; + + cos2 * + (s) " 

- 2 (* + s) S* *«* ' - 1 1 81,129 • ' 2 Wr (' + 5 ) 


e 

n 2 




sin 2(f) 


3 cos 2 6 + 


( 1+ s) 


_ f <f> 


Figure 1 shows the rotating spherical coordinate system to which equations (la-c) refer. 



FIG. 1 Rotating Spherical Coordinate System [9] 


(la) 

(lb) 

(lc) 


where 6 is the in-plane (pitch) angle, <j> is the out-of-plane (roll) angle, and r is the distance from 
the orbital center. The orbital center 0 moves in a circular orbit at constant angular velocity Cl. No 
external forces are applied to the body. Equations (1) are valid for short tether lengths since they were 
derived by linearizing the gravity gradient force. 

In order to obtain a simplified solution of system (1) during the deployment phase, the following 
assumptions are made. The validity of such assumptions (based on considerations drawn from ref. [10]) 
may be verified by an a-posteriori analysis. Then, if <f> is small: 

sin <j) ~ <f> 


The system (1) becomes 


cos <f> — 1 
6 <f> 

n ^ 1 n * 1 


^- = 3cos 2 e 


(2a) 
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(2b) 


8 r 8 3 . . r 

TP + 2 Wrn + 2 + sm * > = - 2 o; 


\ I +2 ^ +(i+3cos2< ' w=0 < 2c > 

During the deployment phase, in general, of course, 0 is relatively large. However, when 9 is small 
(initial instants of motion), eq. (2a) yields an exponential solution for r (i.e., T - = const) and eqs. (2b) 
and (2c) can be easily integrated. Moreover, in this case the solution of (2b) and (2c) is exponentially 
damped, and if such a damping is relatively fast as compared to the time required for the whole 
deployment, eq. (2b) can give a “good” prediction of the steady-state solution for the in plane angle 9. 
Then, eq. (2a) can be integrated again assuming 9 = const = 9 as . If 9 is small we have: 


r = 3fi 2 r 

(3a) 

0 f 0 __ r 

(3b) 

fi2 +2 n^n +3 = ~ 2 H? 

^ + 2 ^ + 4<f> = 0 

(3c) 


These equations are basically analogous to those obtained in refs. [10,11], 


If r(0) = r Q then eq. (3a) yields 

r = r 0 e^ nt 


( 4 ) 


Note that r 0 = \/3Qr 0 . Another solution may be 


t 0 = 0 


( 5 ) 


r = r 0 cos h(V3£li) 


( 6 ) 


- = V3Q tgh (y/3Qt) 

T 

Assuming = \/3, eq. (3b) becomes 


( 7 ) 


^ + 2\/3 ^ + 38 = -2V3 (8) 

For 0(0) = 9 0 = 0, and 0(0) = 9 0 = 0, eq. (8) yields 

^^ = -(l + V3m)e-^ nt (9) 

U3S 

where 9 SS = — = —1.15 rad ~ -66°. 

In conclusion, during the deployment pitch, oscillations damp around the steady-state value of -66°. 
This means that the tether drifts forward with respect to the orbiter motion direction. 

Assuming that the tether reaches the steady-state pitch angle very quickly, the deployment time tj 
can be evaluated by integrating the equation 


= 3cos 2 0 s 


( 10 ) 
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Then we obtain 



where r d is the nondimensional deployment time 
L = 20 km, r 0 = 200 m [8], eq. (11) yields 


1 , L 

— ■= In— 

2v37r cos 66° T o 


( 11 ) 


and L is total length of the deployed tether. Assuming 


r d = 1.04 (12) 

If fi = 1.156 x 10~ 3 rad/s [12], we have 

t d ~ 5655s (13) 

The deployment is accomplished in about 1 orbit. Of course this is not a detailed evaluation of the 
deployment time because, besides the assumptions made in deriving the basic equations, the “drift” 
phase only has been modeled, not including the “straighten” and the “braking” phases in which the 
tether tension (neglected in the previous analysis) is gradually raised. Since 10 minutes after the 
ejection the end mass should be over 200 m away from the orbiter [8], the total deployment time is 
estimated to be approximately 6200 s (r from 0 to 20 km). In ref. [8] the duration of the deployment, 
computed by the BEADSIM program (including straighten and braking effects) is 8100 s. For more 
consideration, see Part III: Conclusions and Recommendations. 

As far as the evaluation of the acceleration is concerned, it has to be remembered that in our model 
the deployment trajectory is practically a straight line inclined 66° forward from the local vertical (see 
also [13]). The only component of the acceleration is then the radial one, 

a r = r= ifl 2 r 0 e v/3nr (14) 


Its maximum value for t = t d is 


r m = = 1.34 x 10 2 m/s 2 (15) 

z 

From results given by ref. [8], in the straighten part of deployment the maximum value of acceleration 
may be estimated to be one order of magnitude less than the value predicted by eq. (15). In conclusion, 
the (radial) acceleration in the deployment phase (at maximum) is estimated to be of the order of 10” 3 p 
(g s= 9.81 m/s 2 ). It has to be noticed that the equation for the out-of-plane angle is of the same kind as 
that for the in-plane angle, however the forcing term is not present; accordingly during the deployment 
phase possible out-of-plane displacements are damped around zero angle. Then, in the present model 
the deployment trajectory may be considered totally contained in the orbital plane. 

IL2 - SWING PHASE 


The swing phase may be basically studied by using the equations of system (la-c). It may be 
assumed that swing phase begins when the deployment is completed. In effect, during the deployment 
phase the velocity dependent forces due to the motion relative to the rotating orbital system (Coriolis 
forces) force the tether forward from the local vertical. When the deployment ends Coriolis effects 
vanish and the tether experiences a wide free libration (swing phase). 


Under the same assumptions made to study the deployment phase, and putting r = 0, eqs. (2b) 
and (2c) yield 


9 

n 2 

± 

fi 2 


— 3 sin 0 cos 0 


— 4 sin 0 cos </> 


(16a) 

(16b) 
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Note that in the aim of decoupling the equations in the second equation, cos 9 — 1. 
Equations (16a-b) can also be derived by writing the equilibrium of the pendulum subjected to the 
gravity gradient field. For initial conditions 

9 — 9 m 9 = 0 

= fim <f> = 0 

the solutions of eqs. (16a) and (16b) are 


^ = ±73 \J sin 2 9 m — sin 2 9 

1 u 

(17a) 

^ = ±2\J sin 2 <j> m - sin 2 <f> 

(17b) 


The oscillation period tg may be computed as 



After resolving the integral, we have 


tjfl 

2n 


73 


ld- 



sin 2 9 m + 




sin 6 9 m + ■ 


tail 1 

2tT y/3 COS Offi 

The duration t 3 of the swing phase is then t s = ^ . In nondimensional form 

_ t s fl _ 1 1 

* 2ir 4 73 cos 9 m 


(18) 


If 9 m = 66° 


Tg = 0.23 


(19) 


Equation (19) means that the swing phase is accomplished in about 1/4 of orbit. In other words 


t s ~ 1250s 


( 20 ) 


As far as the order of magnitude is concerned, results (19) or (20) are in agreement with those of 

refe. [8,12]. Once ^ and $ have been calculated, they may be substituted into complete eqs. (la-c) in 
order to evaluate tension T and accelerations. In the swing phase we have 


T = 3 M0 2 L 


cos 2 9 cos 2 —\ + 

u 




and then 


T = 3MQ. 2 L< cos 2 4> 


cos 2 9 + ^sin 2 9 m — sin 2 9^ ± ^ (sin 2 9 m — sin 2 0) 




( 21 ) 
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From eq. (21) it is possible to evaluate the tension T as a function of 9 m ,<fi m and current 9 and <j> 
angles. 


The most difficult situation occurs when <f> = 0. In this case: 


T = 3 MSl 2 L 


cos 2 9 + (sin 2 9 m — sin 2 0) ± y I (sin 2 9 m — sin 2 9) + 1 sin 2 <f> n 


If 9 = 0 (end of swing phase), 9 m = 66° and (f> m = 30°, we have 

T = 9.66Mfi 2 L (22) 

Assumming fl = 1.156 x 10“ 3 rad/s, M = 50 kg, L = 20 km, (22) yields 

T = 12.9JV. 

Note that [see eq. (21)] for high values of 0 m , and 9 near in magnitude to 0 m , in the retrograde phase 
of a period the tether could become slack. This should not be the case in SEDS. 

The results obtained for angular velocities, swing phase timing, and tension agree well with similar 
results of ref. [14], Transverse acceleration components may be calculated from angular accelerations. 
In spherical coordinates it results [15]: 


a# = L9 cos <j) -f 2 L9cj> sin <f> 

” L9 2 . OJL 

a^-L<t> — sin 2<p 

Of course, in the swing phase, r = 0 and r — L. Then 9 and <j) are respectively 


6 = SI 1 


2 { l + tijti tg(l> - 3 2 sin29 


n 2 

<p = — — sin 20 


3 cos' 


* + H)' 


Putting (24a) into (23a) yields 

glq = n 2 L cos <fi 
and, taking into account (17a) and (17b): 


4( i+2 C) 


tg 0 — - sin 29 


ag = Cl 2 L cos <j) ±4 yj sin 2 4> m — sin 2 (p (l ± 2y/%\j sin 2 9 m — sin 2 <?) tg <j) — ^ si 


sin 29 


In an analogous way it results 

Q 2 L 


H 


■ sin 2<p 


1 + 3 cos 2 9 ± 2'fZyJ sin 2 9 m — sin 2 0 + 6 ^sin 2 9 m — sin 2 9^ 


(23a) 

(23b) 

(24a) 

(24b) 


(25a) 


(25b) 
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Let us examine (25a). Suppose <f> m is small, thus eq. (25a) becomes 

a$ = Sl 2 L ±4 ^ - </)2 (l ± 2v^\/sin^jn _ — sin^^ 4> ~ | sin 20 

The factor - 4> 2 is of order of <f> 2 \ then the relation may be simplified as 

ag = — ^ft 2 Lsin20 (26) 

aQ reaches its maximum value at 9 = — \\ clq max = §fi 2 T. For f) = 1.156 x 10~ 3 rad/s and L = 20 km, 
CLQ ~ 4.0 x 10“ 2 m/s 2 . For 9 = 0,a$ = 0. As far as the expression of a 0 is concerned, note that if 
<f> = 0,(1$ = 0. If <f> ^ 0, assuming the plus sign in eq. (25b) during swing phase, the maximum value 
of a.0 is obtained when 0 = 0. (The tether is vertical, before cutting.) The absolute value is 

a0 = fi 2 L sin 2<f>(2 + V3 sin 0 m 4- 3 sin 2 6 m ) (27) 

If 9 m = 66°, a0 ~ 6fi 2 Lsin2<£. If <f> = 15°, ~ 8.0 x 10” 2 m/s 2 . From these results it is evident that 

maximum transverse out-of-plane acceleration can be about 2 times the maximum value reached at 
9 = — J by the transverse in-plane acceleration. 

However it is interesting to observe that for 4> = 1°, O0 = 5.6 x 10 -3 m/s 2 , that is to say, that even 
in the presence of very small out-of-plane oscillations, the out-of-plane transverse acceleration is “only” 
one order of magnitude less than the maximum in-plane acceleration. 

Part III - Conclusions and Recommendations 

The dynamics of motion of an SEDS has been studied by using a simplified model in which no 
external forces (except the gravity gradient field) are applied on the tethered body and the tether is 
massless. Furthermore, due to the short tether length which is typical of SEDS, equations have been 
derived by linearizing the gravity gradient force. The dynamics of SEDS operation has been modeled 
as a sequence of two phases: the deployment phase and the swing phase. In the first one the velocity 
dependent forces due to tethered end mass motion relative to the rotating orbital system (Coriolis 
forces) force the tether itself forward from the local vertical. A wide in-plane angle is quickly reached 
practically at the beginning of the deployment. When the deployment ends, Coriolis effects vanish 
(f = 0) and the swing phase begins, which is characterized by a wide free libration. The deployment 
phase has been studied assuming that the oscillations of the in-plane (pitch) angle quickly damp and 
a wide steady-state in-plane angle is attained practically instantaneously. Nevertheless, in order to 
evaluate such a steady-state angle, the linearized form of the in-plane motion equation has been solved 
(this equation generally is typical of a nonlinear oscillator with damping and driving terms proportional 
to r/r). (This allowed us to evaluate the constant deployment angle to be about 66° from the vertical, 
forward with respect to the orbiter motion.) Accordingly, the deployment occurs at this constant 
in-plane angle for its whole duration. 

Note that usually the payload is ejected with an initial velocity from the orbiter by means of a 
spring. Due to the particular environmental conditions experienced by the tethered mass, this mass 
does not ” feel” instantaneously the gravity gradient field (allowing the tether to be deployed in a free- 
nul (or low) tension motion), but a certain delay occurs. To evaluate this delay we may use a modified 
form of Froude number 

Fr = V 2 /g'L (28) 

where V and L are reference velocity and length and g' is the gravity gradient field acceleration: 

g' = 3fi 2 r (29) 
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Then 


Pr = V 2 /3 fiV 


(30) 


When the payload is ejected from the orbiter at velocity V , it feels the orbital dynamic 
F = 0 [1]. Since r = Vt, eq. (30) gives 


V3SI 


effects if 
(31) 


If Q — 1.156 x 10“ 3 rad/s, t = 500 s. In conclusion, on timescales longer than about 10 minutes, 
orbital dynamic effects become significant. This result agrees with considerations made in ref. [8]. This 
means that, after payload ejection, in the first 500 s the deployment rate may be assumed constant, 
f = 0.4 m/s (ref. [8]), and a distance r = (0.4 m/s) x 500 s = 200 m is initially attained. Prom this 
instant eq. (10) applies. The deployed length, deployment rate, and acceleration are found to be 
respectively 


T = roe Vzn(cos66°)(t-to) 


(32) 


f = v / 3fi(cos66°)r (33) 

t = 3fi 2 (cos66°) 2 r (34) 


Note that eqs. (32), (33), and (34) overestimate r, f, and r with respect to results given by ref. [8]. 
This is probably due to the fact that, in the presence of control tension and if the elasticity of tether 
material is included into the model, the time constant in eq. (32) is less than y/ 3. 


Maximum estimated acceleration is of order of 10 2 m/s 2 (or 10” 3 g). At the end of the “straighten” 
phase ref. [8] estimates f = 0(10 -3 m/s 2 ). 


Deployment time (from r = r Q ) is 


r d ~ = t=- In— (35) 

2n 27r\/3cos66 t 0 

For = 1.156 rad/s the deployment is accomplished in about 1 orbit. 

Note that the steady-state angle in ref. [8] is estimated to vary about from 65° to 72° at the end of 
the straighten phase. Steady-state value is reached in about (500 + 1000)s. 

The swing phase has been assumed to start from 6 = —66°. In order to study this phase the same 
model of deployment phase was used in which, of course, r = 0. 


Relationships for angular velocity and angular acceleration in in-plane and (if possible) out-of-plane 
motions have been derived. See eqs. (17a-b) and (24a-b). The swing duration is 

t s rt i i 

r 3 = — - = -— ===== ~ 0.23 orbits, or t s = 1200 s (36) 

27r 4 v3 cos 0 m 

A complete equation to evaluate the tension as a function of in-plane and out-plane angles is given by 
(21). When 0 = 0° (end of swing phase) it results 

If = ^ = 0 ° T-11.6N (37) 

If = 30° and <j> = 0° T ~ 12.9 N 

The result of eq. (37) is in excellent agreement with ref. [8]. Prom an analysis of eq. (21) it arises that 
for high values of 0 m , and 9 near to 9 m in magnitude, in the retrograde phase of a period of oscillation, 
the tether could become slack. This should not be the case in SEDS. 
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Transverse in-plane and out-of-plane accelerations are given by eqs. (25a-b); ag is practically 
independent of out-of-plane motion and reaches its maximum value at 9 = : a# max = |ft 2 L ~ 

4 x 10” 2 m/s 2 ;a^ attains its maximum value when 9 — 0°. This maximum value, for 9 m = —66° and 
<f> = 15° is max = 3 fi 2 L. Then, the maximum out-of-plane transverse acceleration could be about 
2 times greater than the maximum value of in-plane transverse acceleration . 

Note that most of the data from ref. [8] are taken from table 8.2.3.2, page 158, and are shown in 
Table 1. 

TABLE 1. BEADSIM Summary Output File for 
First SEDS Experiment [8] 


Time 

NomLng 

ActLng 

TotRng 

XRange 

YRange 

ActdL 

RngRate 

Tension 

0 

1 

1 

1 

0 

-1 

0.40 

0.40 

0.011 

300 

115 

115 

115 

-38 

-109 

0.37 

0.38 

0.011 

600 

231 

231 

231 

-134 

-189 

0.40 

0.40 

0.011 

900 

355 

355 

355 

-262 

-239 

0.42 

0.42 

0.011 

1200 

483 

483 

483 

-402 

-268 

0.43 

0.43 

0.011 

1500 

610 

610 

610 

-539 

-284 

0.42 

0.42 

0.011 

1800 

733 

733 

733 

-669 

-301 

0.41 

0.41 

0.011 

2100 

861 

861 

861 

-794 

-331 

0.44 

0.44 

0.011 

2400 

1005 

1005 

1005 

-929 

-384 

0.53 

0.53 

0.012 

2700 

1185 

1185 

1185 

-1090 

-467 

0.68 

0.69 

0.012 

3000 

1425 

1425 

1425 

-1300 

-583 

0.92 

0.92 

0.013 

3300 

1743 

1743 

1743 

-1581 

-733 

1.21 

1.21 

0.014 

3600 

2157 

2157 

2156 

-1954 

-912 

1.54 

1.55 

0.015 

3900 

2677 

2677 

2676 

-2434 

-1113 

1.91 

1.92 

0.018 

4200 

3312 

3312 

3310 

-3030 

-1330 

2 . 30 

2 .30 

0.021 

4500 

4060 

4060 

4057 

-3746 

-2 557 

2.68 

2.68 

0.025 

4800 

4922 

4922 

4918 

-4580 

-1790 

3.05 

3.05 

0.030 

5100 

5896 

5896 

5888 

-5527 

-2028 

3.43 

3.41 

0.037 

5400 

6975 

6979 

6966 

-6583 

-2279 

3.78 

3.78 

0.045 

5700 

8172 

8172 

8155 

-7745 

-2554 

4.17 

4 .15 

0.056 

6000 

9483 

9483 

9462 

-9016 

-2870 

4.58 

4 . 57 

0.071 

6300 

10930 

10930 

10905 

-10409 

-3253 

5.07 

5.06 

0.096 

6600 

12533 

12533 

12509 

-11937 

-3739 

5.63 

5.65 

0.136 

6900 

14 322 

14 322 

14304 

-13618 

-4375 

6.29 

6.33 

0.212 

7200 

16320 

16320 

16310 

-15446 

-5238 

7 . 00 

7 .03 

0.375 

7500 

18482 

18 4 8 4 

18481 

-17318 

-6451 

7.12 

7.13 

1 . 104 

7800 

19770 

19774 

19774 

-18032 

-8116 

1.60 

1.59 

1.940 

8047 

20000 

20007 

20007 

-16993 

-10560 

0.57 

0.58 

4 . 002 

8100 

20000 

20005 

20005 

-16538 

-11256 

0.23 

0.24 

2 . 873 

8400 

20000 

20016 

20015 

-11899 

-16094 

0.02 

0.01 

7.812 

8700 

20000 

20023 

20023 

-2769 

-19831 

0.01 

0.01 

11.614 


EndSec ApsisAlt PeakTension £ sec BrakeWork Exposure Risk:3d>D 
8700 35.63 km 11.620 new 8699 5788 j 16.57 KmHr 0.182% 
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Concluding Remarks 


A first step to enhance the results obtained analytically is to integrate the same model [eqs. (1)] 
numerically. Then mass of tether, elasticity, and tension control law could be included (together with 
aerodynamic effects). To be investigated: 

1) influence of T/M p (T = tension, M p = payload mass) on libration angle and range rate. 


2) influence of Mt/M p (M t = tether mass) on tether curvature. 

3) possible vibrational modes 
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